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Abstract. We have analysed the full ISO spectrum of the planetary nebula NGC 6302 in order to derive the 
mineralogical composition of the dust in the nebula. We use an optically thin dust model in combination with 
laboratory measurements of cosmic dust analogues. We find two main temperature components at about 100 
and 50 K respectively, with distinctly different dust compositions. The warm component contains an important 
contribution from dust without strong infrared resonances. In particular the presence of small warm amorphous 
silicate grains can be excluded. The detection of weak PAH bands also points to a peculiar chemical composition 
of the dust in this oxygen-rich nebula. The cool dust component contains the bulk of the mass and shows strong 
emission from crystalline silicates, which contain about 10 percent of the mass. In addition, we identify the 92 
/xm band with the mineral calcite, and argue that the 60 /im band contains a contribution from the carbonate 
dolomite. We present the mass absorption coefficients of six different carbonate minerals. The geometry of the 
dust shell around NGC 6302 is studied with mid-infrared images obtained with TIMMI2. We argue that the cool 
dust component is present in a circumstellar dust torus, while the diffuse emission from the warm component 
originates from the lobes. 
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1. Introduction 

All low and intermediate mass stars end their life on the 
Asymptotic Giant Branch (AGB) by ejecting their entire 
H-rich envelope through a slowly expanding, dusty wind. 
After t he AGB, the central otar quickly inoroaooo ito of 



characterized by dense AGB remnants and a very hot lu- 
minous central star. 

NGC 6302 is probably one of the best studied PNe 
with a massive progenitor. A recent determination of the 
mass of the ion ized nebula is about 2 Mq ( Pottasch fc 



fective I toniporaturc to valuco high enough to begin ion 
izing its AGB ejecta: a planetary nebula (PN) is born. 
Depending on the mass and thus luminosity of the star, 
the transition from AGB to PN can go very fast: for the 
most massive objects, in less than 1000 years ionization 
of the nebula begins. These massive objects therefore are 



Bcintema 1999), based o n a distance determination of 1.6 
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kpc ( Gomez et al. 1993 ). Dust mass estimates also indi- 
cate a very massive shell {Aid — 0.02 A/q), using a dis- 
tance of 2.2 kpc ( [Gomez et al. 1989 ). These distance de- 
terminations rely on VLA observations of the expansion 
of the nebula over a 2.75 yr period, assuming an outflow 
velocity of 13 km s^^. In principle, this is a very reli- 
able method, however the epoch over which the nebula is 
observed, is very short. Therefore the increase in size is 
very small and hard to measure, indicating that the er- 
ror bars on these distance determinations are still very 
large. Instead, we adopt a distance of 0.91 kpc, based on 
emission- line photometry of NGC 6302, from which the B 
and V magnitudes, the luminosity and the distance can 
be derived ( Shaw fc Kaler 198S| ). The large number of 
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luminous PNe studied gives confidence in the distance de- 
termination from this statistical method. 

The nebular abundances indicate that NGC 6302 is 
a type I nebula, consistent with a massive progenitor. 
Casassus et al. (2000D estimate that the progenitor mass is 



4-5 Mq. The morphology of the nebula observed at opti- 
cal wavelengths is highly bipolar, pointing to non-spherical 
mass loss on the AGB, resulting in a dusty torus in the 



equatorial region (Lester & Dinerstein 1984). The incli- 
nation angle of the system is ^ 45° with respect to the 
line-of-sight (Bohigas 1994). The temperature of the cen- 
tral star is very high: Casassus et al. (2000) mention a 
temperature of ~ 250 000 K, while [Pottasch et al. (1996| ) 
arrive a temperature of ^ 380 000 K. Although there is 
some uncertainty about the distance and therefore masses 
and luminosities involved, everything points to a rather 
massive and luminous progenitor. 

The Infrared Space Observator y (ISO) spectrum has 



been present e d in se veral papers ( Barlow 1998 ; Molster 



et al. 2001a , 2002b) and is characterized by a wealth 
of narrow solid state features in the 20-70 spectral 
range caused by circumstellar dust in the AGB remnant, 
as well as by strong emission lines from a multitude of 
fine-structure lines originating from the ionized gas in the 
nebula. The dust bands have been identified with crys- 
talline silicates and a number of other com ponents (pK^oike 



et al. |000| ; [Molster et al. 20014 |002^. [Kemper et"al 



(2002b, hereafter referred to as Paper I) have reported on 
the detection of carbonates in the dust shell, based on the 
identification of broad features at ^--^60 and ~92 fim. 

Optical and near-infrared images have already shown 
that the dust distribution around NGC 6302 is rather 
complex (Lester & Dinerstein 1984; ^chwarz et al. 19*9^ ). 
The ISO spectroscopy supports this, because a broad dust 
temperature range is required to explain the shape of the 
spectral energy distribution. In addition the dust compo- 
sition seems to be very complex, as there is evidence for 
a mixed chemistry by the presence of both oxygen-rich 
dust and carbon-rich dust features, the latter in t he form 
Polycy clic Ar omatic Hydrocarbons (PAHs) (see Molster 
et al. [ 2001a ) for a discussion on the origin of this di- 



chotomy) . 

In order to reconstruct the mass loss history of 
NGC 6302, including the geometry and composition of 
the AGB wind, infrared spectroscopy and imaging are 
needed. Unfortunately, the ISO data lack spatial informa- 
tion, limiting the analysis to the bulk dust composition. 
Mid-infrared imaging can reveal the present-day geometry 
of the dust envelope, which puts limits on the mass loss 
history. 

This paper is organized as follows: In Sect. ^ the 
ground-based mid-infrared imaging and spectroscopy 
is presented. The ISO spectroscopy is discussed in 
Sect. H, which describes our laboratory data of carbon- 
ates (Sect. 3.1) and a model fit to the observed spectrum 
(Sect. In Sect. | we propose a possible geometry of 
NGC 6302 and discuss the astronomical relevance of car- 
bonates. Our results are summarized in Sect. ^. 



2. The TIMMI2 observations 

2.1. N- and Q-band imaging 

We have observed NGC 6302 using the TIMMI2 mid- 
infrared imaging spectrograph attached to the 3.6m tele- 
scope at the European Southern Observatory, La Silla, 



Chile. For a description of the instrument see Reimann 
|et al. (2000| ). The observations were carried out on the 
night of June 17/18, 2001. Observing conditions were not 
optimal, with occasional clouds and variable seeing (be- 
tween about 0.6" and 1.4"). Imaging was done in four 
photometric bands, centred at 9.8 (N2), 10.4, 11.9, and 20 
(Q) fj,in. In addition, we obtained two 8.35-11.55 /im long- 
slit spectra. We observed HD 81797 and HD 169916 as 
point sources in order to determine the shape of the point 
spread function, as well as for flux calibration purposes, 
which is only performed in case of spectroscopy. 

To correct for background, chopping and nodding has 
been performed. Apart from the Q-band image, the chop- 
ping and nodding images fell all inside the frame. For the 
various N-band images we used a lens scale corresponding 
to a pixel size of 0.3"/pixel, whereas the Q-band image was 
obtained with a pixel size of 0.2"/pixel. The dimensions 
of the detector are 320x240 pixels of which the central 
300x220 pixels were used for further analysis. The chop 
throw used for the images in the N-band is 30", with a 
nodding offset of 30" perpendicular on the chop direction. 
For the Q-band image, we applied a chop throw of 40", 
and a nodding offset of 40" in the same direction as the 
chop. 

Both the positive and negative images were used for 
the resulting image. The final image has been sharp- 
ened by means of a deconvolution using the point spread 
function (PSF) derived from HD 81797 for the Q-band. 
Unfortunately, for the N2, 10.4 and 11.9 filters no stan- 
dard star was measured in the same night with the same 
settings, and we used observations of standards of another 
night. For the N2 and 10.4 filters no standards were avail- 
able at all. Therefore, we used HD 169916 measured with 
the 11.9 ^m filter for all three filters for the PSF, assum- 
ing that the PSF would be similar in all three N-bands. 
A comparison with other N2-band images (with different 
settings, but the positive and negative images did overlap 
each other) showed that our assumption is not unreason- 
able. 

Fig. shows the final TIMMI2 images at N2, N10.4, 
Nil. 9, and Q band, as well as the "raw" (i.e. before de- 
convolution) N2 image, and an optical VLT image of the 
nebula for comparison ( |ESO 1998[ ). The "raw" N2 im- 
age shows a bright elongated region of approximately 10" 
by 4", and a much fainter region offset to the south and 
west. There is also evidence for faint arcs of emission that 
are reminiscent of the X-shaped nebula seen in the op- 
tical image. In fact, the mid-IR emission has a morphol- 
ogy which is strikingly similar to that seen in the optical 
(see Fig. This is surprising since the mid-IR emission 
originates from (warm) dust, while the optical emission is 
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CN2 raw) 



Fig. 2 J TIMMI2 N2-band image of NGC 6302, overlaid 
with tHe contours of the optical VLT image (Fig. |l]) . 
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Fig. 1.TIMMI2 images of NGC 6302. 
In the upper left corner the VLT -first 
light image of NGC 6302 is shown (|ES0 



1995). The upper middle image is the 
logarithmic non-deconvolved N2 band 
image, which shows the low brightness 
structure somewhat better than the de- 
convolved images. In the upper right 
corner the deconvolved N2 band image 
is shown, in which detailed structure 
in the centre appears. The bottom row 
shows the deconvolved 10.4, 11.9 and Q- 
band images. All images have the same 
orientation and scale. 
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TIMMI2 N-band spectra NGC 6302. The spectra 
tacen 1.2" east and 1.2" west of the brightest (east 
torus) lobe in the image. For comparison the ISO 
n (divided by 5) is given as a dashed line. The 
flux calibration is quite accurate, but the absolute 
bration is only certain up to 40%. The difference 
n of the west and east side of the lobe is also given. 
;s around 9.5 /im are due to telluric lines. 



dominated by ionised gas. The equatorial torus, seen as a 
waist in the optical images, must contain very cold dust 
which does not yet contribute to the flux at 10 /im. We will 
return to this point when we discuss the ISO spectrum. 

The deconvolved N2, N10.4 and 11.9 /im images in 
Fig. ^ are very similar and show a rather "blobby" ap- 
pearance. Note that the faint structure in the background 
is a common artifact of image deconvolution. However, 
the structure in the central area is real, since we see simi- 
lar substructure in different observations deconvolved with 
different point spread functions. Also the 6cm radio map of 



1989). Our images are very similar to the 3.3 /im image 
taken by Casassus et al. (2000), taking into account that 
the orientation of their image is accidentally rotated by 
180° (Casassus, priv. comm.). 

The deconvolved Q band image (Fig. |l|) tends to show 
more spherical emission, although the eastern part of the 
nebula is still dominating the emission. Due to the reduced 
spatial resolution at 20 /im it is not possible to say if the 
emission is equally blobby as seen at 10 /im. 



NGC 6302 shows substructure in this region ( Gomez et al 



F. Kemper et al.: The mineralogy of NGC 6302 



'4 



Fig. 4. The position of the N-band sUt showed on the 
deconvolved N2 image. North is up and east is to the left. 
Both slits were 1.2" x 70" and placed 1.2" east and west 
from the central peak. It is clear that both slit positions 
cover the eastern lobe of the nebula. 



2.2. N-band spectroscopy 

In Fig. 3 we show the spectra taken at the slit positions 
indicated in Fig. 0. We used HD 169916 as a reference 



spectrum, assuming that this star has a F{h') 



flux 



distribution. The comparison with other calibration obser- 
vations showed that the absolute flux calibration is only 
accurate up to about 40%. The relative flux calibration 
is much more accurate. The wavelength calibration pro- 
vided by ESO on the TIMMI2 home page appeared to 
be inaccurate, and we determined a new wavelength cal- 
ibration, based on the [Arm] and [Siv] line, assuming a 
linear pixel to wavelength dependence. The fact that the 
[Navi] line and the PAH features are now at the right 
wavelengths shows that our assumption is valid (using the 
calibration files provided by ESO the 11.3 ^m PAH fea- 
ture was found at 11.7 fim). The comparison with the ISO 
spectrum (dashed line in Fig. 3), gives confidence in the 
new wavelength calibration. 

The spectrum is dominated by the forbidden emission 
fines of [Siv] and [Ann]. However the PAH features at 11.3 
and also 8.6 /im are also found. The noise around 9.5 fim 
is due to telluric lines. 

The east and west slit are both probing the eastern 
lobe of the nebula (Fig. ^ . The difference spectrum is very 
flat, apart from the [Siv] line. The fact that even around 
8.6 and 11.3 /im it is flat, indicates that the distribution 
of the PAHs is quite constant in the bright eastern lobe. 
This is in agreement with the 3.3 fim image obtained by 
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Fig. 5. The ISO 2.4 - 200 fim spectrum of NGC 6302. 



This figure is based on data presented in Molster et al 



(2001a). Clearly visible are the numerous emission lines 



due to fine-structure transitions in the ionized gas. 



3. Analysis of the ISO spectrum 

The ISO SWS and LWS spectra have been presented by 
Barlow (1998| ) and [Molster et al. (2001^ |2002bD and in 
Paper I. We briefly repeat the overall characteristics of 
the spectrum. The SWS and LWS apertures, which have 
sizes in the range 14x20" to 80", cover a large fraction 
of the nebula and therefore provide no detailed spatial 
information. The spectrum, shown in Fig. |^, is dominated 
by a strongly rising spectrum which peaks near 55 /xm, 
and numerous, sometimes rather narrow emission bands. 
Both the continuum and these emission bands are caused 
by the circumstellar dust. In addition, strong, spectrally 
unresolved emission from fine-structure lines are evident, 
which are caused by the ionized nebular gas. These lines 
can be used to determine the chemical abundance of the 



gas and properties of the central star (e.g. Pottasch et al 
1996; Pottasch fc Beintema 1999| ). In this work we focus 
on the analysis of the dust component. 



3.1. Laboratory spectroscopy of carbonates 



Casassus et al. (200C) 



Molster et al. (2001a) provide a useful inventory of 
the solid state emission bands in the ISO spectrum of 
NGC 6302, and we will not repeat their analysis here. 
While many emission bands were successfully allocated 
to (mostly oxygen-rich) dust components, several promi- 
nent emission bands, notably near 29, 48, 60 and 90 ^m 
remained partially unidentified. We have investigated pos- 
sible carriers for the ^--^92 fim band, reported to be present 
in NGC 6302 and NGC 6537 ( [Molster et al. 20014 |2002b| ). 
During the course of our laboratory studies of cosmic dust 
analogues, we noted that the carbonate calcite (CaCOa) 
has a strong resonance at ~92 /im. This prompted us 
to carry out a systematic laboratory study of transmis- 
sion spectra of several carbonate species. A series of car- 
bonates of possible cosmic relevance has been chosen for 
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Table 1. Chemical formula, crystal system and content 
of impurities of the carbonate minerals determined by 
SEM and EDX measurements. The impurities in dolomite, 
magnesite, siderite and ankerite are the carbonates corre- 
sponding to the metals mentioned in the last column. The 
used abbreviations in column 3 describe the crystal system 
of the mineral (trigonal and rhombic). 



Mineral 


Formula 


Crystal 


Impurities 






system 


(wt%) 


dolomite 


CaMg(C03)2 


tri. 


~0.5 Fe, very pure 


calcite 


CaCOa 


tri. 


<0.1 Si, very pure 


aragonite 


CaCOa 


rhom. 


<0.2 Si, very pure 


magnesite 


MgCOa 


tri. 


~ 1.0 Fe; 0.2-0.9 Ca 


siderite 


FeCOs 


tri. 


^ 2.0 Mg; - 5.0 Mn 


ankerite 


CaFe(C03)3 


tri. 


2-3 Mn; 3.5-4.2 Mg 



E 
o 




200 



Fig. 6. Mass absorption coefhcients of natural carbonates 
derived from transmission measurements of the powders 
embedded in KBr and PE. For clarity, the spectra have 
been multiplied by the factor indicated in the plot. 



laboratory investigations, including the natural minerals 
ankerite (origin unknown) , aragonite (Bohemia) , dolomite 
(Eugui, Navarra, Spain), calcite (origin unknown), mag- 
nesite (Steiermark, Austria) and siderite (Schonbrunn, 
Vogtland, Germany). The purity and homogeneity of the 
natural samples are investigated by scanning electron 
microscopy (SEM) and energy dispersive X-ray analysis 
(EDX). The macroscopic phase homogeneity of the car- 
bonate samples was measured by X-ray diffraction. The 
chemical composition and the content of impurities are 
presented in Table |l]. 

The minerals dolomite, calcite, aragonite and magne- 
site turned out to be very pure and homogeneous. The im- 
purities in siderite and ankerite are the corresponding car- 
bonates mentioned in Table ^. These are homogeneously 
distributed within the main mineral component. 



Small carbonate grains with an average size smaller 
than 2 /im have been produced by grinding and sedi- 
mentation in water-free acetone. In order to measure the 
infrared spectrum over a wide spectral range the small 
grains were embedded in transparent material like KBr 
and Polyethylene (PE). Infrared spectroscopy was per- 
formed using a Bruker 113v Fourier transform spectrome- 
ter in the wavelength range between 2-200 /im. The in- 
frared spectral behaviour of the carbonate minerals is 
demonstrated in Fig. ^, and can be found in digitized form 
in the Jena-St. Petersburg Database of Optical Constants 



( penning ct al. 1999| ). 

Magnesite, siderite, ankerite, dolomite, and calcite be- 
long to the carbonates of the calcite group and crystal- 
lize in the trigonal crystal system. Carbonates contain- 
ing cations with larger radii than Ca^"*" crystallize in the 
rhombic crystal system. CaCOa is dimorph and can ex- 
ist in both structures as calcite and aragonite. The typ- 
ical crystalline structure of the carbonates in the calcite 
group is characterized by an alternating sequence of cation 
and anion layers. In both magnesite and calcite, the plane 
carbonate group consists of three oxygen located at the 
corners of an equilateral triangle and a carbon in the cen- 
tre. The bonding between C and O is covalent whereas 
the bonding between the cations and the carbonate anion 
(COa^") is mainly ionic. Each cation is octahedrally co- 
ordinated by six O atoms of the adjacent COs^" ion. The 
crystals of the calcite group with mixed metal composi- 
tion - like dolomite and ankerite - show an alternating 
arrangement of the calcium and magnesium or calcium 
and iron layers. In contrast to the calcite structure the 



C; 



,2+ 



ion in aragonite is coordinated by 9 0-atoms. 



The infrared optical properties are determined by the 
crystalline structure of the minerals. The spectral be- 
haviour of carbonates has been measured several times in 



previous studies (e . g. prlellwege et al. 197G| ; penman 1976 ; 
profino et al. 1997 ; Yamagishi et al. 2001). The measured 



mid-infrared bands arise mainly from the carbonate an- 
ion. The planar structure of this anion and its C3 symme- 
try gives rise to three infrared active modes at about 7, 
11 and 14 /im. The 7 /im band can be attributed to the 
C=0 stretching mode whereas the 11 and 14 /im bands are 
caused by the out-of-plane and in-plane bending modes. 
The low frequency bands beyond 25 /xm are due to trans- 
lation of the metal cations and represent motions perpen- 
dicular and parallel to the plane of the carbonate anions 



(Hellwege et al. 1970) 



3.2. Dust model fit 

3.2.1. Description of the model 

We use the identified solid state components to quantita- 
tively fit the entire ISO spectrum. Inspection of the ISO 
data shows that all dust bands are in emission, suggest- 
ing that the dust shell is optically thin at infrared wave- 
lengths. This substantially simplifies the analysis. Given 
the complex geometry of the nebula, the multiple dust 
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components present in the spectrum, and the still incom- 
plete knowledge of the dust species present, we decided to 
focus our analysis on determining the relative mass con- 
tributions of the known dust species, rather than attempt- 
ing a full 2-dimensional radiative transfer calculation. We 
use the optically thin nature of the dust at mid- and far- 
infrared wavelengths to constrain the dust mass of the 
components using mass absorption coefficients. The thus 
derived dust masses will not depend on the assumed geom- 
etry as long as the emission is optically thin. The problem 
then reduces to deriving the temperature over mass dis- 



In case of a circumstellar dust shell, we may assume 
that the temperature and density distribution follow a 
power law: 



P{r) = Po( — 



(4) 



where r represents the distance to the central star, 
and ro the inner radius of the dust shell. Using dM^^ = 
Airr^ p(r)dr, we can express in terms of the tempera- 
ture distribution in the dust shell, according to 



tribution of the various dust components (e.g. Bouwman 
et al 



OOC). The advantage of using mass absorption co- Md{T{r)) 



efficients in an optically thin dust shell, is that it does not 
require knowledge of the location of the dust. 

We derive the dust mass and temperature distribu- 
tion under the following assumptions: i) the dust shell 
is optically thin at mid- and far-infrared wavelengths; ii) 
all grains have the same size, i.e. 0.1 fxm; iii) all grains 
are spherical; iv) all grains are of homogeneous compo- 
sition; v) the number density of dust grains as a func- 
tion of distance is a simple power law (see Eq. (^, where 
Po — |7ra^pd"-o); vi) the dust species all have the same 
temperature ranges, limited by the temperatures at the in- 
ner radius Tq and at the outer radius Tmax, which are free 
parameters. The result of our analysis is a mass over tem- 
perature distribution Md{T) of all identified dust species. 
In Sect. 



47rporo3 fT{r) 

3-p \n 



(5) 



4.1.1 



we will discuss the validity and the impact 
of these assumptions. 

In the general case of an optically thin dust cloud, 
which specifically applies to an optically thin dust shell, 
the specific intensity emitted by the dust in the cloud can 
be written as 



hiT)^ ( a^B^{T)dV for t<1 
Jv 



(1) 



where is the specific absorptivity. By integrating 
over volume V this can be written as 

I,{T)^N<jMT) or I,{T) ^ AUT)kMT) (2) 

with the absorption cross-section can be described in 
terms of the geometric cross-section and absorption effi- 
ciency, according to cr^ — cFgcomQu- Following for instance 



Bouwman et al. (2000 ), Md{T) is the dust mass at tem- 
perature T. The mass absorption coefficient is given by 
and the number of grains in the volume is represented by 
N. 

When we realize that we can write for the specific fiux 
Fi, = I,y/D^, and under the assumption that the dust 
grains are identical spherical grains with a radius a, we 
find 



1)^ Pd 4 a 



(3) 



with D the distance to the dust cloud and pd the den- 
sity of the material in the dust grain. 



For spherical dust grains the density at the inner radius 
is Po = ^Tra^pdno, with uq the dust grain number density 
at the inner radius, which is different for each identified 
dust component. The density pd corresponds to the den- 
sity of the grain material itself. This can be combined with 
Eqs. (I) and (|) to obtain an expression of the specific flux: 



1 Ana'^ro^no (T{r) 



£|2 



■P 



hB,{T) 



(6) 



where the observable total specific flux Ti, is given by 



I F,{T)dT 

'To 



(7) 



which can be compared directly to the ISO spec- 
troscopy. 

3.2.2. Fit parameters and dust composition 



Following the analysis of [\4olster et al. (2001a| ), we as- 
sumed that enstatite (MgSiOa), forsterite (Mg2Si04) and 
crystalline water ice are present, in addition to amor- 
phous silicates. Diopside ((Ca, Mg)Si03) is believ ed to con- 
tribute to the 60 pm complex ( Koike et al. 2000 ). Diopside 
is in the form of a solid solution between woUastonite 
(CaSiOa) and enstatite (MgSiOs) in the number ratio 
46:54, i.e. Wo46En54. In addition, we considered metal- 
lic iron as a dust component, which is found to be present 
in the dust shells of evolved stars ( Kemper et al. 2002a ), 
and the carbonate species presented in this work. Under 
the assumption that all dust grains are spherical and have 
a radius a of 0.1 /im, we have derived values from lab- 
oratory measurements of amorphous olivine (Dorschnei 
|et al. 1995 ); metallic iron (Henning & Stognienko 1996); 
forsterite and clino-enstatite ( Koike et al. 1999|); diopside 
( Koike et al. 200C| ); water ice (Bertie et al. 1969; Warren 



1984) and carbonates (this work) 

Using Eqs. (^) and (^ the far-infrared spectrum has 
been fitted, under the assumption that all dust species 
are found in the same temperature regimes. Thus, the free 
parameters to be fitted are the number density at the in- 
ner radius uq, which differs for each dust component and 
scales with the mass contained in that component; the 
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Table 2. Empirically derived dust masses by fitting the ISO spectrum of NGC 6302 for the temperature regimes 
30-60 K and 100-118 K. The dust masses are obtained by integrating Eq. (^) over the considered temperature 
range. The absolute accuracy of the dust masses is a factor of three, this is due to the uncertainties in laboratory 
measurements and inaccuracies in the used assumptions. The relative dust masses, i.e. the mass ratios between the 
various identified species, are estimated to have an error of ~50% (see Sect. [1.1. 1| ). The values for the crystalline 
silicate components are lower limits, because we used laboratory measurements that were intrinsically much stronger 
than similar measurements obtained by other groups (Sect. [4.1.1] ). In cases where lower limits are given, we did not 
detect the species indicated. The value indicates the maximum amount that could be present without showing the 
spectral features at a discernible level. 
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iron (or carbon) 
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water ice 
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diopside 
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calcite 
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Fig. 7. The warm and cool dust components required to fit 
the ISO spectrum of NGC 6302. The solid line represents 
the ISO spectrum. The warm (100-118 K) and cool (30- 
60 K) dust components are represented by the labelled 
dotted and dashed lines. 



temperature range defined by Tq and Tmax! and the slopes 
p and q of the density and temperature power laws. For 
the distance D we have used the measured value of 0.91 



kpc ( 3haw fc Kaler 1989 ) and the grain size a is assumed 
to be 0.1 /im. The inner radius ro follows from the val- 
ues of To and q, and is therefore not a fit parameter. If 



p and I are both chosen to be —1/2 the best results are 



obtained. Inspection of the spectrum shows that the tem- 
perature distribution of the dust must be very wide, which 
is incompatible with a simple density distribution of 
the dust. We attempted to fit the spectrum with a single 
range of temperatures but were not successful; apart from 
a cool dust component (30-60 K), peaking near 55 /im, 
we had to introduce a second dust component at 100- 
118 K dominating the spectrum between 10 and 30 /xm 



(see Fig. 0). The results of the fits are shown in Table |[ 
where the mass of each dust component for both temper- 
ature ranges is given. These masses are constrained by no 
which is separately fitted for all identified dust species and 
both temperature ranges, and are obtained by integrating 
Eq. (^ over the considered temperature range. 

The warm dust component could not consist of amor- 
phous olivine, because the strong resonances due to the Si- 
O stretching and 0-Si-O bending modes at 9.7 an 18 /Ltm 
are not observed. Iron is used as the carrier of the warm 
continuum, but as PAHs have been detected in NGC 6302 



( Molstcr et al. 2001a ), amorphous carbon, instead of iron, 
could also explain the continuum emission. A mixed chem- 
istry is not uncommon for evolved stars, it is for example 
seen in IRAS 09425-6040 ( [Molstcr et al. 2001b[) . We con- 
clude that the composition of the warm dust is different 
from the composition of the cool dust. 



The exact composition of the cold amorphous dust is 
unknown. Amorphous olivine is used to produce the re- 
quired continuum emission, but other abundant species 
with a smooth absorption spectrum at far-infrared wave- 
lengths, such as pyroxene and iron, could be contributing 
as well. The mass fractions of forsterite and enstatite are 
lower limits. We used the laboratory data leading to the 
best spectral match with the observations, i.e. those of 



Koike et al. (1999D 



The previously unidentified 92 iim band (Molstcr et al 



2001a) can be fitted very accurately with calcite (CaCOa, 
trigonal), see Fig. |^. The strong resonances of calcite at 
shorter wavelengths (see Fig. ^ are not discernible in the 
ISO spectrum, indicating that it is not a significant com- 
ponent in the warm dust. Aragonite (CaCOa, rhombic) 
also has a resonance at ~92 /im, and may therefore con- 
tribute to the observed feature in NGC 6302. However it 
is not very likely to exist in a Planetary Nebula. Although 
it has the same chemical composition as calcite, its lattice 
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Fig. 8. Model fit to the ISO spectrum 
of NGC 6302. The thin line represents 
the observed spectrum, the fit, consist- 
ing of the components given in Table ^ 
is indicated with the thick line. In the 
lower part of the plot, the contribu- 
tions of the cold dust components (30- 
60 K) are indicated, except the con- 
tribution of amorphous olivine, which 
leads to a smooth overall contribution. 
The inset shows the observed spectrum 
(thin line) and the model fit (thick line) 
from 60-120 fim, where the contribu- 
tions of all species, except calcite, are 
subtracted from both the observations 
and the model. Thus, it only shows the 
observed and modelled calcite 92 /im 
feature. Figure adopted from Paper I. 



Jones 1998) and therefore, we concentrated our analysis 



200 



150 



^ 100 - 




on dolomite. 

The results of our model fit are shown in Fig. |^ and 
Table Since the dust shell is optically thin, the relative 
mass fractions can be derived from the model fit. The ab- 
solute dust mass of each mineralogical component scales 
with the distance, and if 910 pc is adopted for the dis- 



tance (Shaw & Kaler 1989), we arrive at the dust masses 
as indicated in Table We derived that the depletion of 
calcium into the identified calcium bearing species - cal- 
cite, dolomite and diopside - is about 30% (Paper I), if 
we assume a dust/gas mass ratio of 1/100 and adopt a 
calcium abundance of [Ca/H] = 2.2 • 10"^ ( [Snow 1984| ). 
From the fit to the 60 fj,m complex the fraction of water 
contained in the solid phase can be calculated. Adopting 
a mass ratio of 1 • 10""^ between H2O and H2 (Gonzalez 



Fig. 9. 1 Contributions to the 60 ura complex in NGC 6302. 



The thin line represents the ISO spectrum of NGC 6302 
from 40-80 /im with all contributions subtracted, except 
the contributions of diopside, dolomite and water ice. The 
thick lines represent the contributions of diopside, water 
ice and dolomite to F^, indicated with abbreviations of 
their mineral names. The uppermost thick line is the sum 
of the contributions of diopside, dolomite and water ice. 



structure is different, and it can only exist under relatively 
high pressure and temperature. 

The 60 fim feature can be accurately fitted using ab- 
sorption efficiencies of dolomite, diopside and water ice 
(see Fig. |9|). Ankerite (FeCOs) has a resonance at approx- 
imately the same wavelength (see Fig. ^ and could in 
principle contribute to the 60 fim complex. However, so- 
lar system carbonates are almost exclusively in the form 



Alfonso fc Cernicharo 1999] ) we find that ^ 10% of the 
water is contained in water ice. This is consistent with 



the results of Hoogzaad et al. (2002 ). The presence of OH 
molecules ( Payne et al. 1988 ) indicates that water should 
also be present in the vapour phase. 

4. Discussion 

4.1. Geometry and composition of the circumstellar 
dust shell 



In Sect. 3.2.2 we derived that the dust in the circumstel- 



of dolomite and calcite (Johnson & Prinz 1993; Brearley 



lar environment of NGC 6302 consists of two temperature 
components. The cold component (30-60 K) has a mass 
of 0.050 Mq and contains silicates, water ice and carbon- 
ates, indicative of a oxygen-rich chemistry. The warm com- 
ponent (100-118 K) contains only 1.3 • lO^'' Mq, which 
is mainly caused by a featureless species like carbon or 
iron. About 5% of the warm component is in the form 
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warm dust 
PAHs 
scattered light 



cool dust 
OH maser 



Fig. 10. Proposed geometry of NGC 6302. A dense torus 
of cool dust is present around the central star (light grey) 
under an inclination of 45° with respect to the line-of- 
sight. More diffuse material in the polar regions, is indi- 
cated in dark grey. The torus obscures the polar region 
on the western side. A part of the torus is located behind 
the eastern lobe. The emission from the warm dust and 
PAHs originates from the lobes, as well as scattered light 
from the central star. The torus contains the cool dust and 
the OH maser. For clarity the approximate position of the 
central star is indicated, however in reality it is obscured 
by the dusty torus and can only be studied in scattered 
light from the polar regions. 



of silicates. The large difference in chemical composition 
between the warm and the cold component and the lack 
of continuity in the temperature distribution suggest that 
the two components have a different spatial distribution. 
This then allows a mixed chemistry; regions of C-rich dust 
could exist in a PN with a strong silicate spectroscopic sig- 
nature. The presence of PAHs in NGC 6302 (first reported 
as UIR bands by Roche & Aitkcn 1986) supports the idea 
of a mixed chemistry, both the warm, diffuse component 
and the PAHs could reside in the same regions of the neb- 
ula. 

The large mass of the cold component suggests that it 
is present in a circumstellar torus, which leaves the polar 
regions for the warm carbon-rich component. Optical im- 
ages of NGC 6302 show that the eastern lobe is brighter 
than the western lobe (see Fig. |l| and also Schwarz et al 



uggesting that the western lobe is partly obscured 



1992(D 

by thel circumstellar torus. This leads to a geomotr}^ as 
indicated in Fig. ^ The dusty torus around the central 
star is seen under an inclination angle ~ 45°, and partially 
obscures the western lobe. In return, the eastern part of 
the torus is located behind the eastern outflow. The dusty 
torus contains a high column density of dust towards the 



central star, and therefore, most of the dust is effectively 
shielded from heating by the central star and remains rel- 
atively cool. The dust in the polar regions however, is 
directly radiated by the central star and thus has a higher 
temperature. 

Our observations, as well as previously published work, 
is in agreement with the proposed geometry. The torus 
is optically thick in the UV, optical and near-infrared. 
Therefore, the optical images only show the scattered 
light from the central star in the polar regions. We have 
convolved the warm and cool dust spectra modelled in 
Sect. |3.2.2| with the TIMMI2 Nil. 9- and Q-band transmis- 



sion profile and corrected for the atmospheric transmission 
in order to determine the contribution of the dust com- 
ponents in both bands. We find that in the N11.9-band, 
only 0.2% of the flux comes from the cold component, 
and consequently 99.8% originates from the warm dust 
component. For the Q-band these numbers are 29.0% and 
71.0% respectively. Since the peak in intensity in the N- 
band images and the optical image coincide (see Fig. ^ 
we conclude that indeed the warm dust is located in the 
polar regions. The Q-band imaging reveals a different pic- 
ture; one third of the flux originates from the cool dust 
component, and two third from the warm dust. Both the 
warm dust in the eastern lobe, and the cool dust in the 
western part of the torus contribute to the observed im- 
age, leading to a more or less spherical structure in this 
wavelength band (Fig. [|). 

N-band spectroscopy o f the eastern lobe sh ows PAH 
emission. From Fig. 6 of ( Casassus et al. 200C , figure is 
rotated by 180°) it is clear that the PAH emission inten- 
sity peak coincides with the intensity peak in the optical 
and N-band, indicating that the PAH emission also arises 
from the lobes, where UV radiation from the central star 
can penetrate. The UV radiation field also explains the 
ionization observed in the lobes (Fig. 3). [Persi et al. (1999 ) 
also notice that the emission of warm dust grains and for- 
bidden line transitions emerges from the same region. The 
6cm data presented by Gomez et al. (198^ ) traces the lo- 
cation Hii region. In their Fig. 3 the region confined by 
the torus is clearly visible. 

On the other hand, the optically thick torus around the 
central star effectively shields the molecular material in 
the torus from the high energy radiation. The infrared ra- 
diation field in the torus will pump OH molecules present 
and cause the formation of OH masers. Indeed, OH maser 
emission is observed in NGC 6302, on the west side of the 
object (Payne et al. 1988), indicating that it arises from 



part of the torus located at the front side of the object. 
This is consistent with the fact that predominantly the 



blue-shifted part of the line profile is seen (Payne et al 



1988) 



4.1.1. Validity of model assumptions 

The geometry described here is consistent with the as- 
sumptions we have used to fit the ISO spectrum, and with 
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the TIMMI2 images. Here we will discuss the validity of 
the assumptions used in this model. 

i) We assumed that the dust is optically thin at large 
wavelengths. The TIMMI2 images indeed suggest that this 
is true: the torus is still visible as a dark lane in the N- 
band images but starts to disappear when NGC 6302 is 
observed in the Q-band. At longer wavelengths, where the 
spectral energy distribution (SED) peaks and the different 
solid state species are identified, the dust torus will be 
optically thin. 

A simplified calculation verifies this: The dust is 
thought to be present in a torus, with a scale height h 
from the equatorial plane and an outer radius -Rout- The 
inner radius is small compared to the outer radius and 



absorption efficiency Qi,/a. For amorphous olivine, metal- 
lic iron and water ice we used optical constants, and de- 
rived absorption efficiency for grains with a size of 0.1 
/im. Only the mass determination of these three species 
is probably somewhat affected by our assumption on the 
grain size. 

iii) We assumed that all grains are spherical. Again 
an assumption that is probably not very realistic. The 
dominant dust component (amorphous olivine) must be in 
the form of non-spherical particles to explain the infrared- 
spectrum of AGB stars, and the same is true for metallic 



iron ( Kemper ct al. 2002a ) . The similarity in the spectral 
appearance of the solid state features detecte d in PNe 
and their d ir ect progenitors, AGB stars, (e.g. Sylvester 



can be ignored. The torus is seen under an inclination an- ct al. 1999; Molster et al. 2002b ), indicates that the 



gle of 45° ( [Bohigas 1994[ ) . The length of the line-of-sight 
through the disk is thus / 2h^/2. Assuming that the 
density is constant throughout the disk, we find for the 
infrared optical depth along this line-of-sight: 



TiR ^ Qn CTgoom 2hV2 

where 

Md 1 



^^grain 2/l7ri?o 



(8) 



(9) 



By adding all cold dust components listed in Table ||, 
the dust contained in the torus is determined to be Afd — 
0.05 Mq. The mass of a grain follows from the used grain 
size of 0.1 /im, and the average density of silicate grains. 
Using cgs units, we find that the torus is optically thin in 
the mid- and far-infrared in case 



Ront^ 



< 1.1 • 10-3^cm-2 



(10) 



For a typical value of Qir in the mid- and far-infrared 
of 0.0028 we find that for i?out > 5.0 • lO^^ cm, or 3.4 • 10^ 
AU, the dusty torus is optically thin at these wavelengths. 
This size is not unreasonable for a circumstellar torus, as 
it is comparable to the sizes of the dust shell shells around 
AGB stars, the PN progenitors. In addition, we can argue 
that the density distribution will not be flat, but probably 
behaves more like a power law. Then, most lines-of-sight 
through the torus will be optically thin even in case of a 
smaller torus. Therefore, we are confident that the torus 



dust grains in the two types of objects probably have 
similar properties. Therefore, we may conclude that the 
grains in NGC 6302 are very likely not homogeneous 
spheres, as these particles result in a shift in peak po- 
sition with respect to non-spherical particles and inhomo- 
geneous spheres. However, in our model the only purpose 
of the grain shape is to calculate the total mass of some 
of the dust components. The effect that grain shape has 
on light scattering is not taken into account, as our model 
uses mass absorption coefficients. 

iv) We assumed that all grains arc of homogeneous 
composition. Using mass absorption coefficients to calcu- 
late the emerging spectrum, we simply derived the total 
mass of a certain mineral found in the line of sight. It is im- 
possible to discriminate between homogeneous grains, and 
grains consisting of two or more mineral species, therefore 
we cannot exclude the presence of grains of mixed compo- 
sition. For instance, it is well known that volatiles, such 
as water vapour, cond ense in layered structures on grains 
( Jura fc Morris 1985 ). We stress that the derived total 
masses are independent from the assumed homogeneity of 
grains. 

v) We assumed that the density and temperature gra- 
dients are simple power laws of the distance to the central 
star. The complex structure of NGC 6302 indicates that 
this is probably far besides the truth, but since the optical 
depth is low, all infrared radiation emitted by the grains 
is received. Therefore not the density and temperature as 
a function of distance to the central star, but the mass- 
temperature relation of the different dust species becomes 
important, in a similar fashion as described by Bouwman 



is inde ;d optically thin in the mid- and far-infrared. 



ii) We assumed that all grains have a size of 0.1 /im. 
In reality, the grains formed in the circumstellar environ- 
ment will cover a wide range of sizes. However, we believe 
that we actually get a good estimate of the actual dust 
masses, because we use the mass absorption coefficients, 
and is thus independent from grain size. Using the mass 
absorption coefficients will lead to a direct determination 
of the dust mass, without calculating the mass of individ- 
ual grains. The calcite and dolomite data presented in this 



et al. (2000 ). Our simulation of this relation by a power 



law density and temperature gradient seems to be in agree- 
ment with the observed spectrum. 

vi) Finally we assumed that all dust species are found 
in the same temperature range. However, if the different 
dust components are not in thermal contact with each 
other, i.e. if they are present in separate - but co-spatial 
- grain populations, the difference in optical properties, 
notably in the UV and visual, causes a difference in the 



work a "e mass absorption coefficients , the data used for 2001). The determination of the temperature range is 



temperature profile (Bouwman et al. 200C; Kemper et al 



diopside, enstatite and forsterite are measured in terms of based on the overall shape of the SED, and we forced 
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the crystalline components - which are the carriers of the 
narrow features superposed on the broad continuum - to 
have the same temperature profile as the amorphous com- 
ponent. Since the dust emission strongly depends on the 
temperature of the grains, this will probably be our main 
source of error in our mass estimates. However, dropping 
the constraint that all dust species are found in the same 
temperature range, leads to a large increase in free param- 
eters for our simple model. 

Considering all the effects discussed here, we estimate 
that our absolute mass determinations are accurate within 
a factor of two, if the distance toward NGC 6302 is 910 
pc. If we take into account the large uncertainty in the 
distance determinations, the error in the absolute dust 
mass determination further increases. The uncertainty in 
the various laboratory measurements introduces an addi- 
tional uncertainty of a factor of two, leading to a total 
error of a factor of ~3. This inaccuracy in the laboratory 
measurements can be inferred from comparison between 
pu blished data. The c rystalline silicate features published 
by Koike et al. (199£ ) are intrinsically much stronger (up 
to 5 times), than the same features in measurements by 



other groups (Servoin & Piriou 1973; Steyer 1974; Jager 
et al. l |)98| ). Consequently, the crystallinity, defined by the 
mass of enstatite and forsterite as a fraction of the total 
dust mass, lies within 5% - 25%, and the values given in 
Table ^ are thus just lower limits. The error in the rela- 
tive dust mass fractions, i.e. the mass ratios between the 
identified minerals, is much better, as they all depend in 
the same way on our model assumptions. We therefore 
estimate the error in the relative dust mass fractions of 
the order of 50%. The error in the absolute and relative 
flux levels is small compared to the inaccuracies discussed 
here, and is therefore not taken into account. 

In order to further approach reality, full radiative 
transfer calculations should be performed. Using spatial 
information derived from imaging and measurements of 
outflow velocities it is possible to reconstruct the geome- 
try of the nebula. With 2-dimensional full radiative trans- 
fer calculations over a large wavelength range, from UV 
to far-infrared, it is possible to resolve the temperature 
and density distribution throughout the PN, for the dif- 
ferent dust components. Also, the effect of grain sizes, 
grain shapes and mincralogical homogeneity then becomes 
apparent in the spectroscopy. This will lead to a much 
more accurate description of the properties of the dust 
in the nebula. However, we chose not to do these calcu- 
lations for the following reasons: First, optical constants 
are not available for some of the identified dust compo- 
nents. Second, our knowledge about the geometry of the 
nebula is still too limited to justify the huge effort of 2- 
dimensional radiative transfer calculations. These calcu- 
lations introduce new free parameters, and it will be very 
hard to constrain the model parameters in an unambigu- 
ous way with the limited information available. Spatially 
resolved spectroscopy of the mid- and far-infrared regions 
will certainly improve the image we have of NGC 6302. 



Finally, to obtain reliable mass estimates the distance to- 
ward NGC 6302 should be more accurately known. 



4.2. The astronomical relevance of carbonates 

In paper I and this work, we present the first extrasolar 
detection of carbonates. The detection is based on far- 
infrared features of carbonate at 92 /Ltm and dolomite, 
which contributes to the 60 /Ltm complex. Due to the 
low abundances, the strong resonances present at short 
wavelengths (see Fig. ^ do not stand out in the spec- 
tra, as many other species contribute significantly to the 
opacity at these wavelengths. This is not the case at the 
far-infrared wavelengths, therefore it is possible to detect 
even very small amounts of carbonates, beyond reasonable 
doubt. 

There have been a number of searches for carbonate 
features in the mid-infrared. The previously unidentified 
infrared (UIR) feature at 11.3 /im seen in pl anetary neb- 



ula NGC 7027 was attributed to carbonates ( Gillett et al 
|1973| ; [Bregman &: Rank 1975| ). However, this identification 



became unlikely with the non-detection of the 6.8 /im car- 
bonate feature in NGC 7027 ( [Russell ct al. 1977bD , and 
the lack of characteristic features due to carbonates in the 



22-35 //m region of the spectrum ( McCarthy et al. 1978 ). 
The UIR feature at 11.3 /im in HD 44179, also known as 
the Red Rectangle, was attributed to the carbonate mag- 



nesite (MgCOa) (Bregman 1977). However, Cohen et al 
(1986| ) showed PAHs are the most likely carrier of the UIR 



feature at 11.3 /im seen in various astrophysical environ- 
ments. Thus, the sole detection of the 11.3 /xm feature is 
no longer considered evidence for the presence of carbon- 
ates. 

Interesting environments to search for carbonates are 
star forming regions and young stars. Carbonates form 
easily on the surface of planets when liquid water is 
present, but othe r formation mechanis ms cannot be ex- 
cluded (Paper I). Russell et al. (1977a ) have determined 
an upper limit to the carbonate/silicate mass ratio in the 
interstellar medium toward Orion. They find that the this 
mass ratio is at most 0.05, based on the detection limit of 
the 6.8 /im feature. The 6.8 /im feature is seen in absorp- 
tion toward embedded protostars. The resemblance of the 
spectral appearance of the 6.8 /im feature in W33A and 
in Interplanetary Dust Particles (IDPs) has lead to the 
conclusion that they have the same carrier, most likely 
carbonates ( ^andford fc Walker 1985 ; Tomeoka fc Buseck 
1986), although Tielens et al. (1984) argue that it is due to 
hydrocarbons. Cox (198S| ) observes the 6.8 /im feature to- 
wards protostar AFGL 961 and claims that it is partially 
due to carbonates along with a contribution of hydrocar- 
bons. Hints of the 11.3 and 13.4 /im carbonate features 
are also claimed to be present. The similarity between the 
6.8 /im features in IDPs and in interstellar lines-of-sight 



is reported by Quirico et al. (200C) as well, presumably 
coming from the same carrier, although carbonates are 
rejected as a carrier. From the high resolution ISO SWS 
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spectroscopy it becomes clear that the spectral shape of 
the 6.8 /xm absorption feature observed towards protostars 
cannot be explained by carbonates, but a satisfactory al- 
ternative is still lacking ( Keane et al. 200l] ). Recently, the 
~92 /zm carbonate feature has been detected towards pro- 
tostar NGC 1333-IRAS 4 (CeccareUi et al, in prep.). 

Although extraterrestrial carbonates are quite often 
found in meteorites and IDPs, they are actually not easy 
to detect by means of astronomical observations. The only 
claim still standing at the moment is the detection of fea- 
tures at 26.5 and 31 /im due to dolomite and calcite respec- 



tively in the ISO SWS spectrum of Mars (LcUouch et al 
2000()7]H owever, the far-infrared features of dolomite and 
calcite at ~60 and ^^92 /im are not seen on Mars (Burgdorf 
et al. 2 |000|) . 

As pointed out in Paper I, until now carbonates are 
believed to be formed through aqueous alteration. In this 
formation process, carbon dioxide (CO2) dissolves in liq- 
uid water and forms carbonate ions (COa^"). If cations 
like Ca^"*" or Mg^+ are found in the solution as well, car- 
bonates can be formed as a lake sediment when satura- 
tion is reached. The presence of carbonates is seen as ev- 
idence of planet formation, as an atmosphere containing 
carbon dioxide is required, as well as liquid water on the 
surface of the planet. The connection with planet forma- 
tion is the justification for the search for carbonates in 
the circumstellar environment of young stars. Carbonates 
are ubiquitous in our own solar system, as is shown from 
the composition studies of meteorites and interplanetary 
dust particles. Carbonate-containing meteorites are con- 
sidered tr acers of the formatio n history of th e solar sys- 
tem (e.g. Endress et al. 199(: ; Bischoff 1998 ). However, 
the detection of carbonates in PNe suggests that alterna- 
tive formation mechanisms exist and that the presence of 
carbonates no longer provides direct evidence for planet 
formation (Paper I). 



5. Summary 

We have determined the composition and distribution of 
the dust in the planetary nebula NGC 6302. We found that 
a warm (100-118 K) and a cool (30-60 K) dust component 
are present. The cool dust component is located in a cir- 
cumstellar torus of which the inner part effectively shields 
the UV/optical and near-infrared radition from the cen- 
tral star. The torus contains amorphous olivine, forsterite, 
clino-enstatite, water ice, diopside, dolomite and calcite. 
Outside the solar system, the carbonates dolomite and cal- 
cite are only seen in NGC 6302 and NGC 6537 (see also 
Paper I), and the formation mechanism of carbonates in 
these environments is not yet understood. Diopside is so 



far only foun d in NGC 6302 ([Koike et al. 2000| ) and in two 
OH/IR stars (Dcmyk et al. 200C). Water ice, forsterite, en- 
statite and amorphous olivine are very common in the dust 
shells of evolved stars (e.g. [Molster et al. 20024 |b|Jc|). The 
circumstellar torus also contains the OH maser reported 
by [Payne et al. (1988| ). 



The optical depth in the polar regions is smaller, and 
therefore, these regions can be observed in the optical in 
scattered light from the central star. These regions con- 
tain the warm dust component, which thermally emits in 
the mid-infrared. The N-band images are dominated by 
thermal emission from the warm dust. From the ISO spec- 
troscopy, constraints on the dust composition could be de- 
rived: predominantly metallic iron (or another featureless 
dust component), amorphous olivine, forsterite and clino- 
enstatite are detected. For water ice, diopside, calcite and 
dolomite we have derived upper limits for the mass frac- 
tion in the warm component. Only 0.3% of the total dust 
mass is contained in the warm component, which has a 
different mineralogical composition than the cool compo- 
nent. UV radiation can penetrate the polar regions and 
thus we are able to expl ain the observed PAH emission, 
which is also reported by Casassus et al. (2000 ). 

In addition, the mass absorption coefficients of carbon- 
ate minerals are presented in this work. We studied calcite, 
dolomite, ankerite, aragonite siderite and magnesite. 
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